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This paper describes a phase principle jor measuring the location or 
the spectral shape of a discrete radio source. The phase principle is relatively 
simple to implement and leads to a measurement of location or spectral 
shape which is insensitive to receiver gain fluctuations. For measuring 
the location of a weak, discrete radio source, the theoretical accuracy is 
slightly better than the theoretical accuracy resulting from the Ryle inter- 
ferometer. For measuring the spectral shape of a weak, discrete radio 
source, the theoretical accuracy is slightly better than the theoretical ac- 
curacy resulting from either the Ryle interferometer or the Dicke radiometer. 
Furthermore, the implementation of the phase principle doesn't require 
input switching. Also, the calibration curve associated with the phase 
principle is independent of changes in the average receiver gains of the 
two receivers. 

I. INTRODUCTION 

In many branches of science and technology observations of a dis- 
crete radio source provide fundamental knowledge. In the field of 
radar the illuminated target serves as the discrete radio source. In the 
field of space exploration the radio transmitter on-board the space 
vehicle serves as the discrete radio source. In the field of radio astro- 
nomy the "radio star" serves as the discrete radio source. 

The "radio star" is a remarkable example of a discrete radio source. 
In the past twenty years radio astronomers have discovered that na- 
ture provided many discrete radio sources or radio stars at certain 
locations in the sky. What are the locations of these radio stars? What 
is the power spectrum of the observed radiation from a particular 
radio star? Answers to such questions are of fundamental importance 
in the field of radio astronomy. 1 
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In order to measure relatively small values of radiated power from 
a discrete radio source, one must compete with the inevitable back- 
ground noise and the inevitable radio receiver noise. It is well known 
that one requires a method of measurement which is relatively in- 
sensitive to receiver gain fluctuations. The papers by Dicke 2 and 
Ryle 3 discuss this important point in more detail. In fact, the present 
day method for measuring relatively small values of radiated power 
from a discrete radio source makes use of some form of the Dicke 2 
radiometer or the Ryle 3 interferometer. 

The purpose of this paper is to describe a phase principle for meas- 
uring the location or the spectral shape of a discrete radio source. We 
shall see that the phase principle is relatively simple to implement 
and leads to a measurement of location or spectral shape which is 
insensitive to receiver gain fluctuations. For measuring the location 
or spectral shape of a weak, discrete radio source, we shall see that the 
theoretical accuracy associated with the phase principle is slightly 
better than the theoretical accuracy associated with the Dicke radi- 
ometer or the Ryle interferometer. We shall also see that for measur- 
ing the location or spectral shape of a weak, discrete radio source 
using only phase information, the accuracy associated with the phase 
principle is essentially equal to the accuracy associated with the maxi- 
mum likelihood principle. 

II. MEASUREMENTS BASED ON THE PHASE PRINCIPLE 

2.1 Implementation 

Fig. 1 illustrates a simplified implementation of the phase principle 
for measuring the location or spectral shape of a discrete radio source. 
S(t), Niit), and N 2 (t) represent zero mean, independent, narrow-band, 
stationary Gaussian processes. iVi(f) and N 2 (t) each represent the sum 
of background noise plus receiver noise. N t (t) and N 2 (t) are assumed 
to have equal variances. The spacing, d, between the two antennas is 
many wavelengths in order that Ni(t) and N 2 (t) can be considered as 
independent processes. S(t — At) and S(t + At) are due to the presence 
of a discrete radio source located at a small angle with respect to 
boresight. We assume that the receivers preserve the phase difference 
between the antenna excitations. 

S(t — At) + Ni(t) and S(t + At) + N 2 (t) represent the outputs 
of the two receivers. i\i represents the ith independent sample of the 
phase difference between S(t - At) + N^t) and S(t + At) + N 2 (t). 
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Fig. 1 — Simplified implementation of a phase principle for measuring location 
or spectral shape of a discrete radio source. When source is located at d sin = 
\/4, t)i = 0. When source is located at 6 = 0, receivers scan in frequency and ijt 
traces out a measure of spectral shape. 

Hi is taken to be in the primary interval (— t, it). After n such samples 
the output rf is given by 



t 1 f. 



n — 



COS 77, , 



(I) 



where 



n = Bt 

B = IF bandwidth 

t = observation time. 

We shall assume that n is relatively large like n ^ 10 4 , since we are 
primarily interested in observing relatively weak, discrete radio sources. 
Fig. 2 illustrates a relatively simple method for generating t? 1 from 
the inputs S(t - At) + JV,(0 and S(t + At) + N 2 (t). B l , ta , and 0, 
represent the envelope, IF angular frequency, and phase angle, respec- 



S(t-At) + N,(t)=R,cos[(i>ot+0,] 



S(L+At)+N 2 (t) = R 2 cos[o> t + fl 2 ] 



*V2 cos[w o t+0i] 



cos(0 2 -0,) 




^VFcos[<i) t+ff 2 ] 



Fig. 2 — A method for generating r;t from the inputs S{L—\t)-\-N,(l) and 
S(t-\-At)-\-N?(t). The band-pass limiters remove all amplitude information. 
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tively, of the narrow-band Gaussian process S(t — At) -f- Ni(t). Sim- 
ilarly, R 2 , co , and 6 2 represent the envelope, IF angular frequency, 
and phase angle respectively of the narrow-band Gaussian process 
S(t + At) + N 2 (t). The band-pass limiters shown in Fig. 2 are well- 
known devices for removing all amplitude information and preserving 
the phase information as is indicated in Fig. 2. See Davenport and 
Root 4 for a discussion of the band-pass limiter. 

As indicated in Fig. 2, if one takes the product of a/2 cos (o) t + 0,) 
and \/2 cos (<a t + 6 2 ) and passes the result through a suitable low- 
pass filter the result is cos (0 2 — #0 = cos tj(£). By taking the average 
of this result, one can generate r^. 

This method of generating q* can also be used to help implement 
the phase principle described in Ref. 5 in order to detect the presence 
of a discrete radio source located at 6 = 0. 

Fig. 2 indicates clearly that ^ is independent of receiver gain fluctua- 
tions or changes in the average receiver gain of each receiver. Also, 
the two receivers need not have the same average gain. Thus, an 
unusually long observation time r is advantageous when using the 
phase principle. 

Notice that if the band-pass limiters in Fig. 2 are shorted out, we 
have the well-known correlator configuration. ' 7,80 

We shall now show that a measurement of r? leads to a measurement 
of location or spectral shape of the discrete radio source. 

2.2 Statistical Properties of 77,- 

In order to simplify the analysis, we shall always assume that the 
discrete radio source is at a small angle with respect to boresight. 
To begin, we shall state some known statistical properties of the angle r?,- . 

Equation (34) of Ref. 10 gives the probability density p 2 (i?) of each 
independent sample 17 ,• as 

!>*(V) ■- -^- (1 " &T l \k ^n- ft + *p + VT-W,~\ , 

where 
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Var = Variance 

7/0/2 = — d sin — o) At 

A 

7T . -1 T 

The Fourier series development of (2) follows from Middleton's" 
equation (9.33) 

f\ 1 x 1 ^ /T 2 (n/2 + 1) „(n n , , A , v 

(3) 
where 2^1 is the Gaussian hypergeometric function 

»/ «. . v __ -, , a& , a(a + 1)0(0 + 1) x 2 
Jfafiiyix) = 1+—X + t(t + 1} 2!+ ••• 

and r is the gamma function. 

The expectations E cos 77, and E cos 2 77,- follow from (3) 

E cos t/,- = j 2 F X {\, \;2;l 2 ) cos ti« (4) 

E cos 2 * - 5 + f 2^(1, 1; 3; I 2 ) cos 2^ . (5) 

We shall see that the phase principle for measuring location or spectral 
shape of a discrete radio source is based on (4). Equation (4) should 
be compared with (4) of Ryle, 3 the equation which characterizes the 
output of a Ryle interferometer. Both equations are proportional 
to cos 779 . 

2.3 Measurement of Location 

Let us first consider the problem of measuring the location of a 
discrete radio source whose true location is some small positive angle 0. 
From (4) we see that E cos 77,- = when rj g = r/2 or d sin = X/4. 
This suggests that we observe T7 + and conclude that d sin = X/4 
when 77 + = 0. How accurately can we form an estimate of in this 
manner? 

For 7} e near x/2 let the estimate yje of 779 be determined from the 
linear equation 

v = j 2^(1, t; 2; 0(| - *)• (6) 
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Thus, 



Var fj g = (2tt)-( - 1 Var 6 = 
or, in a more suitable form, 



Var rj 



(f^(M;2;J 2 )J (7) 



n(^) 2 Var = 4 (1 + a) 2 (tf cos 2 vd^ih h 2; 0, (8) 

where E cos 2 ??,- is given by (5) with t] B = tt/2. Equation (8) char- 
acterizes the theoretical accuracy associated with the phase principle 
for measuring the location of a discrete radio source and is plotted in 
Fig. 3. 

2.4 Measurement of Spectral Shape 

Now let us consider the problem of measuring the spectral shape 
of a discrete radio source located at 6 = 0. We shall assume that the 
variances of the background noises and receiver noises are invariant 
over the frequency region of interest. Under these conditions the esti- 
mate d of the signal-to-noise power ratio "a" can serve as an estimate 
of spectral shape by using the well-known frequency scan technique 

0.24 



dsiN0=-£- 




0.4 



Fig. 3 — Theoretical accuracies of Ryle interferometer, phase principle, and 
maximum likelihood principle (using only phase information) for measuring 
the angular location of a discrete radio source. 
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indicated in Fig. 1. How accurately can we form an estimate a of "a" 
in this manner? 

Equation (4) with = or tji = defines "a" as an implicit func- 
tion of E cos 77, which we shall indicate by 

a = H(E cos 7i,). (9) 

Equation (9) suggests that we form an estimate a of "a" from the 
equation 

d = H( v ). (10) 

A plot of (10) is presented in Fig. 4. This figure can be considered as a 
theoretical calibration curve. Notice that the theoretical calibration 
curve is independent of changes in the average receiver gain of each 
receiver. This is indeed unusual. One measures ip and reports the cor- 
responding value of d. Assuming that Var 2Vi and Var N 2 are invariant 
with frequency over the frequency range of interest, d will then trace 
out the spectral shape of the discrete radio source as the receivers scan 
in frequency. We shall now characterize the accuracy of the estimate &. 
For large n, the only case of interest in this paper, CrameYs 12 work 
shows that the estimate d is characterized, approximately, by a Gaus- 
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Fig. 4 — The theoretical calibration curve associated with the phase principle 
for measuring spectral shape of a discrete radio source. The receivers scan in 
frequency and a traces out the spectral shape. 
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sian probability density having the following expectation and variance: 

Ed = H(E cos Vi ) + 0(0 = a (11) 

Var d = H\ Var j + 0(n l ) = H\ Var n \ (12) 

where 

dd [ dE cos Vi I ~| -1 

coa„ ~ L da loJ 

= *£^t-£ [| ,*•,(§, f ; 3; + Aft. *; 2; J 2 )]" 

This last result follows by differentiating (4) with respect to "a", 
setting tj e = 0, and then taking its reciprocal. Equation (11) implies 
that the estimate d is essentially unbiased for large n. 
From (12) and (1) we have 

n Var & = H\(E cos 2 tj, - E 2 cos t,,), (13) 

where E cos ij,- and E cos 2 ^ are given by (4) and (5) with rj e = 0. 

Equation (13) characterizes the theoretical accuracy associated with 

the phase principle for measuring the signal-to-noise power ratio "a" 

or the spectral shape of the discrete radio source and is plotted in Fig. 5. 

III. MEASUREMENTS BASED ON THE RYLE INTERFEROMETER OR DICKE 
RADIOMETER 

3.1 Measurement of Location 

When 6 is small and t? 9 = t/2 or d sin 6 = X/4, the theoretical ac- 
curacy associated with the Ryle interferometer for measuring the 
location of the discrete radio source was derived by Manasse. 18 In 
our notation Manasse's 13 (60) becomes 

n(yf Var B = (2tt)- 2 (1 + a) 2 . (14) 

Equation (14) characterizes the theoretical accuracy associated with 
the Ryle interferometer for measuring the location of the discrete radio 
source and is plotted in Fig. 3. 

3.2 Measurement of Spectral Shape 

One can measure the spectral shape of the discrete radio source 
located at 9 = by using the Ryle 3 interferometer or the Dicke 2 radi- 
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Fig. 5 — Theoretical accuracies of Ryle interferometer, Dicke radiometer, phase 
principle, and maximum likelihood principle (using only phase information) for 
measuring the spectral shape of a discrete radio source. Both the Ryle inter- 
ferometer and the Dicke radiometer require some amplitude information. 

ometer. In fact, these methods are at present the accepted methods. 
We shall go on to characterize the accuracy associated with these 
methods of measuring spectral shape. 

In order to simplify the notation of this section, let P. = Var S(t) 
and P N = Var JV,(0 = Var JV 2 (Z). Then, 



p. , , A 

a — -=r- and a = :=-• 

* N "n 



(15) 



P, denotes the unbiased estimate of P. and P N is regarded as a param- 
eter. 

The Ryle 3 interferometer utilizes a phase reversing switch to produce, 
periodically, the following inputs to a square-law detector 



25(Q + A r ,(A + N 2 (t) 
2 



(16) 



or 



2 



(17) 
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Thus, using some of Rice's 14 results, the mean values at the output 
of the square-law detectors are, periodically, 

2P > + Pn (18) 



or 

?f. (19) 

The difference in the outputs of the square-law detector is taken as 
the unbiased estimate P. . Thus, EP. = P. , and Ed = a. 

Also, using some of Rice's 14 results, the variances at the output of 
the square-law detector are, periodically, 

[2P, + P N ] 2 , 20) 

4(n/2) l 

or 

[P*l 



4(»/2) (21) 

The factor n/2 appears because the output is in either position only 
one half of the time. Since the difference in the outputs of the square- 
law detector is taken as the unbiased estimate P. , the variance of 
P, is given by the sum of expressions (20) and (21): 

YttP , . pp. + pf + w«r (22) 

2n 

or 

n Var & = nP~ N 2 Var P. = 2a 2 + 2a + 1 . (23) 

Notice that the Ryle interferometer can be considered as a Dicke 2 radi- 
ometer switching between the two inputs given by expressions (16) 
and (17). Thus, (23) characterizes the accuracy of both the Ryle 
interferometer and the Dicke radiometer for measuring the spectral 
shape of the discrete radio source. Equation (23) is plotted in Fig. 5. 

IV. MEASUREMENTS BASED ON THE MAXIMUM LIKELIHOOD PRINCIPLE 
USING ONLY PHASE INFORMATION 

4.1 Measurement of Location 

If one uses the maximum likelihood principle 12 to process a large 
number n of independent samples of the phase difference 77, in order 
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to estimate the location of the discrete radio source when rj e 
or d sin 6 = X /4, one finds 



n \\) Varfl= w a Var * 



= tt/2 



(24) 



I 



(2tt) _ U_, _p a Vo d7?„/ _ T/2 



1 /l dp. 



f/77 



where p 2 is given by (2) and the integrand of the integral in equation 
(24) is to be evaluated at t? 9 = t/2. For various values of "a", the 
definite integral appearing in (24) was evaluated numerically by using 
a digital computer and Simpson's rule. The resulting curve is plotted 
in Fig. 3. Incidentally, this curve applies for all values of 779 . 
As a -» we find that (8) and (24) both yield 



lim 



' (ad 



Yar 6 



= — = 0.02053. 



(25) 



Thus, as a — > 0, the phase principle and the maximum likelihood 
principle using only phase information are essentially equivalent. 

4.2 Measurement of Spectral Shape 

If one uses the maximum likelihood principle 12 to process a large 
number n of independent samples of the phase difference 77, in order to 
estimate the signal-to-noise power ratio "a" or the spectral shape of a 
discrete radio source when = 0, one finds 



n Yar d = 



Jo Pi Lda J 



n' V 



(26) 



where p 2 is given by equation (2) with rje = 0. For various values of 
"a", the definite integral appearing in (26) was evaluated numerically 
by using a digital computer and Simpson's rule. The resulting curve 
is plotted in Fig. 5. This curve also applies for all values of t} g . 
As a -> we find that (13) and (26) both yield 



lim [n Yar d] = -^ = 0.81057. 

,.-0 IT 



(27) 



Thus, as a — > the phase principle and the maximum likelihood principle 
using only phase information are essentially equivalent. 
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V. COMPARISONS OF THEORETICAL ACCURACIES 

5.1 Measurements oj Location 

When using the Ryle interferometer to measure the location of a 
weak, discrete radio source located at d sin 6 = X/4, we have, from (14), 

lim \n(^\ Var d] = ^ = 0.02533. (28) 

Whereas, when using the phase principle or the maximum likelihood 
principle to measure the location, we have, from (25) , 

lim [?i(y) 2 Var §j - 4 = 0.02053. (29) 

Thus, the phase principle and the maximum likelihood principle are 
essentially equivalent, and they are both slightly more accurate than 
the Ryle interferometer. 

See Fig. 3 for a comparison of the theoretical accuracies at other 
values of "a". 

5.2 Measurements of Spectral Shape 

When using the Ryle interferometer or the Dicke radiometer to 
measure the signal-to-noise power ratio "a" or the spectral shape of a 
weak, discrete radio source located at 6 = 0, we have, from (23) , 

lim [nVard] = 1. (30) 

o-«0 

Whereas, when using the phase principle or the maximum likelihood 
principle to measure the signal-to-noise power ratio "a" or the spectral 
shape, we have, from (27) , 

lim [n Var d] = ^ = 0.81057. (31) 

Again, the phase principle and the maximum likelihood principle are 
essentially equivalent, and they are both slightly more accurate than 
either the Ryle interferometer or the Dicke radiometer. 

See Fig. 5 for a comparison of the theoretical accuracies at other 
values of "a". 

For values of "a" away from zero, Fig. 5 shows that the Ryle inter- 
ferometer or Dicke radiometer are more accurate than the maximum 
likelihood principle using only phase information. Thus, one must 
conclude that the Ryle interferometer or the Dicke radiometer require 
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some amplitude information. Consequently, their accuracy is subject 
to deterioration by gain variations. 

Notice that (29) divided by (28) equals 8/ir 2 , and (31) divided by 
(30) also equals 8/tt~. Thus, for measuring the location or the spec- 
tral shape of a weak, discrete radio source, Var and Var a associated 
with the phase principle are lower, by the same factor 8/tt, than the 
corresponding variances associated with the Ryle interferometer. 

VI. CONCLUSIONS 

For measuring the location or the spectral shape of a discrete radio 
source, the phase principle leads to a measurement which is insensitive 
to receiver gain fluctuations. 

For measuring the location or the spectral shape of a weak, discrete 
radio source, the accuracy associated with the phase principle is 
slightly better than the accuracy associated with the Ryle interferom- 
eter or the Dicke radiometer. Also, the accuracy associated with the 
phase principle is essentially equal to the accuracy associated with 
the maximum likelihood principle using only phase information. 

The phase principle is relatively simple to implement, and the im- 
plementation doesn't require input switching. 

The calibration curve associated with the phase principle is inde- 
pendent of changes in the average receiver gain of each receiver. The 
two receivers need not have the same average gain. 

An unusually long observation time is advantageous when using the 
phase principle. 

For measuring spectral shape, both the Ryle interferometer and the 
Dicke radiometer require some amplitude information. Consequently, 
their accuracy is subject to deterioration by gain variations. 
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